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Abstract 

Motivated by the study of matter waves in Bose-Einstein condensates and coupled 
nonlinear optical systems, we study a system of two coupled nonlinear Schrodinger 
equations with inhomogeneous parameters, including a linear coupling. For that 
system we prove the existence of two different kinds of homoclinic solutions to the 
origin describing solitary waves of physical relevance. We use a Krasnoselskii fixed 
point theorem together with a suitable compactness criterion. 

Key words: Nonlinear Schrodinger systems, solitary waves, fixed point theorems 
in cones 



1 Introduction 



Nonlinear Schrodinger (NLS) equations appear in a great array of contexts 
[53P9] as for example in semiconductor electronics p^1[30] . optics in nonlin- 
ear media [27], photonics [25, plasmas [T7j, fundament at ion of quantum me- 
chanics [H] , dynamics of accelerators [19] , mean- field theory of Bose-Einstein 
condensates [Hj or in biomolecule dynamics [THj. In some of these fields and 
many others, the NLS equation appears as an asymptotic limit for a slowly 
varying dispersive wave envelope propagating in a nonlinear medium [17] . 

In this paper we will be interested in systems of two coupled NLS equations 
in which there is a linear coupling between both components. 
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A first example which arises in the study of a spinor Bose-Einstein conden- 
sate composed of two hyperfine states (for instance in many experiments, the 
\F = l,mf = —1) and \F = 2,m/ = 1) states of ^''Rb atoms are used) and 
coupled by an optical or r.f. field. These systems have received a lot of atten- 
tion both experimentally |32ll33|35|36|31|llj and theoretically (see e.g. Refs. 
[55|20|21 II 16|26|l9|34y3 71145] and references therein) since they represent the sim- 
plest mixture of different ultracold quantum degenerate gases and were studied 
immediately after the historial achievement of Bose-Einstein condensation in 
1995 tl4j. 

To simplify the treatment yet preserving the spatial aspects we assume the 
condensate to be magnetically tightly confined along two of the transverse 
directions to a single effective dimension [39] . In the mean-field approximation 
the system is described by the Gross-Pitaevskii equations 

= (Li + Un\H^ + UulH'') + ^^2, (la) 

= (L2 + U22\H' + U2l\^Pl\') ^2 + A^i, (lb) 

where Lj = —d'^/dx^ + Vj with j = 1,2. Equation ([1]) is written in di- 
mensionless form: the spatial coordinates and time are measured in units 
of i = \J hj {2muj) and l/tu, respectively, while the energies and frequencies 
are measured in units of hoj and tu, uj being the trap frequency in the (x, z)- 
plane. The dimensionless nonlinear coefficients, for the quasi-one-dimensional 
condensate, are given by Uij = Uji = 2Naij/i where the scattering 

lengths for binary collisions, N is the total number of atoms, and i is the 
oscillator length in the y-direction. The normalization of the wave function 
ip=co\{ipi,ip2) is then / ip'^ip(Pr = 1. In many practical situations for spinor 
condensates we have Uu ~ f/12 = f/21 — f^22 and it is customary to consider 
all of the scattering lengths to be equal. 

An specific situation which arises also in applications to BEG happens when 
two different condensates do not coexist but are spatially separated and weakly 
coupled. In that situation, described for example in Refs. [12||29] and references 
therein the model is equivalent to Eq. [T]but with different nonlinearities given 
by U12 = U21 = (i.e. there is no nonlinear coupling between both equations). 

Although the dynamics of coupled Bose-Einstein condensates has attracted 
a lot of attention in the last years. Equations ([T]) also arise in other physi- 
cal contexts, most notably in Nonlinear Optical models, where they describe 
coupled optical fibers where the functions Ui and U2 describe the light field 
within each waveguide [56]. Also other physical situations in Nonlinear Optics 
are described by this type of model equations [27j . 
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In this paper we will consider the physically relevant question of the existence 
of solitary wave solutions to models such as the one given by Eqs. ([1]) but 
in more general scenarios where all the parameters (i.e. the potential, the 
nonlinear coefficients and the coupling coefficients) are functions depending 
on the spatial variables. 

The study of the properties of localized solutions and propagating waves 
in Bose-Einstein condensates with spatially inhomogeneous interactions has 
been a field of an enormous level of activity in Physics in the last few years 
[15|5U|4ip!|22|54|HS]|71l38lll8|4U|[5] motivated by novel experimental ways to con- 
trol experimentally the interactions (i.e. optical manipulation of the Feschbach 
resonances [5T]). 

In this paper we will consider theoretically the problem of the existence of 
solitary wave solutions in two-component condensates when all of the physical 
parameters are spatially dependent, i.e. the potentials Vj{x), the nonlinear 
coefficients Uij{x) and the coupling parameter A(x). The latter posibility has 
not been explored yet, neither experimentally nor theoretically. First because 
the analyses of spatially dependent interactions have been focused on the case 
of a single component |l3|50|41p!|22ll54ll46|f7|38|18|40|[8] . Secondly because the 
coupling coefficient A is usually taken to be spatially homogeneous. However it 
is very simple to add an spatial dependence to this parameter in experiments 
since A is related to the coupling field and in the optical case this means just 
using a spatially dependent laser, i.e. a beam with a prescribed profile on the 
scale of the BEC. 

Thus we will consider the physically relevant problem of studying solitary 
wave solutions to the problem 
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' dt = 


dx"^ 



+ a{x) 
+ d{x) 



ipi - h{x)ip2 - c(a;)F(V'i, ^2), 
^2 - e{x)ipi - f{x)G{ipi, 'tp2)- 



(2a) 
(2b) 



where ipj are the complex wavefunctions defined on all M which for solitary 
waves, i.e. localized solutions, must also decay at infinity and we consider 
general nonlinearities F and G satisfying certain conditions to be made precise 
later. In order to look for stationary solutions we eliminate time as usual by 
defining 

iljj{x,t) = Uj{x) exp^iujt). (3) 

Since all physically relevant nonlinearities are homogeneous in time and re- 
defining the potentials a{x),d{x) to incorporate the constant factors propor- 
tional to UjUj we get the final equations to be studied in detail in this paper 
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— Ui{x) + a{x)ui{x) — b{x)u2{x) = c{x)F{ui, U2)ui, (4a) 
—U2{x) + d{x)u2{x) — e{x)ui{x) = f{x)H{ui, U2)u2- (4b) 

i.e. a set of two coupled stationary linearly coupled nonlinear Schrodinger 
equations with spatially inhomogeneous coefficients. We will keep in mind 
the previous discussion from which it is clear that the two most interesting 
situations from the point of view of applications are 

F{ui,U2) = ul, H{ui,U2) = ul (5a) 

for weakly coupled Bose-Einstein condensates or optical fibers and 

F{ui,U2) = H{ui,U2) = uj + ul, (5b) 

for multicomponent spinor Bose-Einstein condensates coupled by an optical 
field [55j- The purpose of this paper is to study the existence of localized 
positive solutions u = (^1,^2) of Eqs. (jlj) satisfying 

lim Mj(x) = 0, lim u^(x) = 0, i = 1,2, (6) 

>oo 1x1— ►oo 

where a, b, c, d,e, f E L°°(M), in Eqs. (HI), are non-negative almost everywhere 
and F{ui, U2) and G(mi, M2) are continuous functions. Such solution u will have 
finite energy, that mathematically means that it belongs to if^(M) x H^{W). 

Between the different solutions of Eq. (j4]) those with a more direct physical 
interest are the so called ground states, which are positive real solutions of 
that system. Although we will be mostly interested on those solutions we will 
also consider the existence of other families of solutions which are relevant for 
applications. The existence of positive solutions for different types of vector 
Nonlinear Schodinger Equations has attracted a lot of interest in recent years 
from the mathematical point of view. For instance in Refs. [3]|^ll5]|6] the ex- 
istence of positive solutions was proven using either critical point theory or 
variational approaches. 

In that context in this paper we complement previous studies with the analysis 
of the richer system given by Eqs. (HI), which describes physical situations 
beyond those studied previously. To do so we will also develop a new theoretical 
approach to the problem based on the use of a fixed point theorem due to 
Krasnoselskii for completely continuous operators defined in cones of a Banach 
space together with a suitable study of the Green's function for the linear 
part of the problem. This method has been employed successfully for scalar 
problems on the real line and in some problems on bounded domains 
[T2] . Using this technique, we will be able to prove the existence of positive 
solutions of Eqs. (jlj) under the conditions IQ. We think that, beyond its 
own interest, this topological approach appears to be complementary to the 
variational approach. 
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The rest of the paper is organized as follows. In Section [2] some preliminary 
results are collected. Section [3] contains the main result about the existence of 
a homoclinic orbit to the origin (positive solution). In Section S] we prove that, 
if in addition to the hypotheses of the main result of Section 12] we assume a set 
of hj^otheses on the functions a, b, c, d, e, f when there exists a second wave 
which is odd. Finally, section [5] contains a result about the study of branches 
of solutions dependent on a parameter. 

We use the notation M+ = (0,+oo), M+ = [0,+oo).For a given a E L°°(M), 
the essential infimum is denoted as a*. The support of a given function a is 
denoted by Supp(a). The limit value of a given function u in +oo (or — oo) is 
written simply as m(+oo) (or u{—oo)). 



2 Preliminaries 

The proof of the main results is based on a well-know fixed point theorem 
in cones for a completely continuous operator defined on a Banach space, 
due to Krasnoselskii [28]. We recall the statement of this result below, after 
introducing the definition of a cone (see, for example, |23j). 

Definition 1 Let X be a Banach space and P be a closed, nonempty subset 
of X . P is a cone if 

(1) \x + fxy E P Vx, y E P and VA, /i G M, 

(2) X, —X E P implies x = 0. 

We also recall that a given operator is completely continuous if the image of 
a bounded set is relatively compact. 

Theorem 1 Let X be a Banach space, and let P d X be a cone in X . Assume 
Qi, Q2 OjTC open subsets of X with G fii, f2i C ^2 and let T : Pn{VL2\Vti) P 
be a completely continuous operator such that one of the following conditions 
is satisfied 

(1) \\Tu\\ ^ ifu E PndVti, and \\Tu\\ ^ if u E Pndn2. 

(2) \\Tu\\ ^ \\u\\, ifu E Pndni, and \\Tu\\ ^ \\u\\, if u E PndVL2. 

Then, T has at least one fixed point in P (1 {Q2\^i)- 

This result has been extensively employed in the study of nonlinear equations 
[52|57] and also in the study of boundary value nonlinear systems p!2|25] . 
However, for problems defined in non-compact intervals such as ours, there 
is the difficulty that the Ascoli-Arzela theorem is not sufficient for proving 
the complete continuity of the operator. We shall employ the following com- 
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pactness criterion (reminiscent of a result by M. Zima [57]) to show that the 
operator is completely continuous. 

Proposition 1 Let Q C BC{M.). Let us assume that the functions -u G are 
equicontinuous in each compact interval of M and that for all u & Q we have 



\u{x)\<^{x), VxgM (7) 



where ^ G BC(R) verifies 



Then, Q is relatively compact. 



lim ^(x) = 0. 

|a;|— >+oo 



Proof: Given {un}n a sequence of functions of Q, we have to prove that there 
exists a partial sequence which is uniformly convergent to a certain u. Note 
that the elements of Q are uniformly bounded by ||^||oo and equicontinuous on 
compact intervals by hypothesis, therefore the Ascoli-Arzela theorem provides 
partial sequence (call it again {un}n) which is uniformly convergent to a certain 
u on compact intervals. Of course, u satisfies also ([7]). Now, we have to prove 
that 

We > 0, 3no s.t. n > uq =^ \\un — m||oo < £■ 

By using ([8]), fix A; > such that max < e/A. On the other hand, 

xm\]-k,ki I 

by using the uniform convergence on compact intervals, there exists no such 
that max \un[x) — u[x)\ < ej^ for all n > no. Then, 

ll^n — m||oo< max \un(x) — u(x)\ + max \un(x) — u(x)\ < 

x£[-k,k] xm\]-k,k[ 

<- + 2 max \f(x)\<e, 
2 xm\]-k,k[ ^' 

and the proof is finished. ■ 

In order to apply Theorem [H we need some information about the properties 
of the Green's function. Let a G L°°(R), a* > 0. For the homogeneous problem 



-(f)" + a{x)4 
M-oo) = 0, 



= 

foo) 



0, 



(9a) 
(9b) 



the associated Green's function is given by 
Gi{x,s) = 



0l(x)02(s), — OO < X < S < +00 
4>l{s)(p2{x), — OO < S < X < +00 



(10) 



where 0i,02 are solutions such that 0i(— oo) = 0, 02(+oo) = 0. Moreover, 
01,02 can be chosen as positive increasing and positive decreasing functions 
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respectively. For a given h{x) G ^^"^(M), the function u{x) = G{x, s)h{s)ds is 
the unique solution of the equation —0" + a{x)(f) = h{x) in the Sobolev space 
if^(M). In particular, u and u' vanishes at ±cx). See for instance 



Note that 0i,02 intersect in a unique point Xq. So, we can define a function 
Pi E BC{R) by 

;i/^.w..<.o, 

Pi{x) = { (11) 
l/0l(x), X > Xq. 



The following result was proven in 

Proposition 2 The following properties for the Green's function defined by 
M) hold. 



(PI) Gi{x, s) > for every {x, s) G M x M. 

(P2) Gilx, s) < Gi{s, s) for every {x, s) eRxR. 

(P3) Given a non-empty compact subset P C M, we define 

mi(P) = min(0i(inf P),02(supP)). (12) 

Then, 

Gi{x,s)>m{P)p{s)Gi{s,s) for all (x, s) G P x R (13) 

Obviously, the linear operator L2[u\ = —u" + d{x)u with d^ > defines a 
second Green's function G2{x, s) with analogous properties involving properly 
defined P2{x), m2(P). 



3 Existence of positive bound states 



From now on, we will assume that M = Supp(6) USupp(c) USupp(e) USupp(/) 
is a non-empty compact set. In order to apply Theorem [T], we take the Banach 
space X = BG{M.) x BG{R) with the norm ||u|| = maxj=i_2|^j|; ioT u = 
(mi,M2) G X. Let us define the cone 

P = \u = (ui,U2) E X : Ui(x),U2ix) > Vx, mm Ui(x) > rniPnlluAl > 

{ x&M J 

(14) 

where pi = miMPi{x), pi = miMP2{x) and the constants mi = mi(M) and 
m2 = m2(M) are defined by property (P3). Note that the compactness of M 
implies that pi > and pg > 0. Moreover, as from (P3) it is easy to see that 
rriipli < 1, for i = 1,2. Thus, one can easily verify that P is a cone in X. 
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Let T : P — s> X be a map with components (Ti, T2) defined by 

Ti(u)(x) = / Gi{x,s) (b{s)u2{s) + c{s)F{ui,U2)ui{s))ds 

= / Gi{x,s) {b{s)u2{s) + c{s)F{ui,U2)u2is))ds. (15a) 

JM 

T2{u){x) = / G2{x,s) {e{s)ui{s) + f{s)H{ui,U2)u2{s))ds 

JR 

= / G2{x,s) {e{s)ui{s) + f{s)H{ui,U2)u2{s))ds. (15b) 

JM 

A fixed point of T is a solution of system (jll) wliicli belong to H^(R) x if^(M), 
and therefore verifies the boundary conditions ([6]). 

Lemma 1 Let us assume that 

F{ui, U2), H{ui, U2) ^ for every ui, U2 ^ 0. 

Then, T{P)CP. 

Proof: Take u = (^1,^2) G -P- The property (PI) of the Green's function 
together with 6*, c* > implies Tiu(x) > for all x. Let us call Xm the point 
where minx^MTiu is attained. Then, for all x G M, 

Ti{u){xm) =Ti{ui,U2){xm) = / Gi (x^ , s) (6(s)m2 (s) + c(s)P(mi , ^2)^1 (s) ) 

>mi / p(s)Gi(s, s) (6(s)m2(s) + c(s)P(mi, ■U2)mi(s)) rfs 

>mipl / ^1(3;, s) (6(s)m2(s) + c(s)P(mi, U2)ni(s)) (is 

= miPoTi(ui,M2)(a;) = miPoTi(u)(a;) (16) 

where we have used (P2) and (P3). In a similar way, we can prove that 
T2{u){xm) > '"^2Po^2(u)(a:;). This completes the proof. ■ 

Lemma 2 T : P —>■ P is continuous and completely continuous. 

Proof: The continuity is trivial. Let us prove that the components of T are 
completely continuous. Let Q C P he a, bounded set, with C > a uniform 
bound for its elements. The functions of Ti{Q) are equicontinuous on each 
compact interval (in fact, the derivative is bounded in compacts). On the 
other hand, for any u G 

\Ti(u)(x)\<C [ GAx, s)b(s)ds + C max F(ui,U2) [ GUx, s)c(s)ds =: ^(x). 

jR l|u||<C JR 
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Since the supports of b, c are compact, ^ G BC(R) fl ^^^(M), therefore Ti{Q) is 
relatively compact by Proposition [H The proof for T2 is analogous. ■ 

The following one is the main result in this section. 

Theorem 2 Let us assume the following hypotheses, 

(a) M is a non-empty compact set. 

(Hi) M2), H{ui, U2) ^ /or etJer?/ tii, U2 ^ 0. 

(^ivj There exist tq > anc? 7, A; > sfic/i that given < r < tq, 

'U2), H{ui, U2) < kr'^, for all ||u|| < r 

(^fj There exist Rq > and 6,K > such that given R> Rq, 

F{ui,U2), H{ui,U2) > KR\ if ui e [miPo-R, -R] or U2 G [m2plR,R]. 

(vi) /^j s)b{s)ds, /^j ^2(3^, s)e(s)(is < 1 /or a// x G M. 

Then, there exists a non-trivial solution u G X of the system 

Proof: We define the open sets Qi and Q2 as the open balls in X centered in 
the origin and with radius r and R, respectively, to be fixed later. Let us take 
u ePndni. Thus 

||Tu|| = max^eR ( / Gi{x, s) (c(s)F(mi, U2)ui{s) + b{s)u2{s)) ds, 

/ G2{x, s) {f{s)H{ui, U2)u2{s) + e{s)ui{s)) ds ) 

JM J 



\Jm JM J 

+ rmax( / Gi{x, s)b{s)ds, / G2{x, s)e{s)ds] < r = ||u|| (17) 

^-■eiR \Jm JM J 

for r sufficiently small, where we have used the hypothesis {iv) and {vi). 

On the other hand, let us take u G P fl dVL2. Thus, at least a component 
of u, say Ml, satisfies UmiHoo — R (the case ||m2||oo — -R is similar). Then, 
rriiplR ^ Ui{x) ^ R for all x G M. Therefore, by using hypothesis (v) we get 

||Tu|| = max (|Tu|)j > max / Gi{x, s) {c{s)F{ui,U2)ui{s) + b{s)u2{s)) ds 

i=l,2,x£R xdM JM 

> KR^+^ max / GAx, s)c(s)ds >R = ||u|| (18) 

xeM J M 

for R sufficiently large. For the case ||'U2||oo = R) everything works in the same 
way. 
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Now, Theorem [T] guarantees that T has a fixed point u G P fl (f22\l^i). Thus 
< < -R so it is a non-trivial solution. ■ 

Although the hypothesis (vi) appears to be quite technical, its meaning can be 
clarified by realizing that the function /^^ Gi{x, s)b{s)ds (resp. /^^ G2{x, s)e{s)ds) 
can be interpreted as the unique solution of the equation —u" + a{x)u = b{x) 
(resp. —u" + d{x)u = e{x)) with boundary conditions u{—oo) = = n(+oo). 
Having this idea in mind, we can formulate the following corollary. 

Corollary 1 The later result still holds if (vi) is replaced by 
(v% J max , , < 1, max——— < 1 

xm a{x) a:GK d{x) 

Proof: Defining u{x) = Jj^j Gi{x, s)b{s)ds, we need to prove that ||m||oo < 1- 
Take xq such that u{xo) = ||m||oo, then u is convex in a neighborhood of xq, so 
from the equation —u" + a{x)u = b{x) we easily get 

«(x„) < ^1^4 < 1. 

a{xo) 

In the same way, J^^ G2{x, s)e{s)ds < 1 and Theorem [2] applies. ■ 



4 Odd solitary waves in the equation with symmetric coefficients 



The aim of this section is to extend some ideas presented in [52l Section 4] 
for a scalar NLS equation to the vectorial case. In the previous section, we 
have obtained positive solutions of the problem under consideration, which is 
the case considered in the all related papers known to the authors. In this 
section, we prove the existence of a new kind of solution when the coefficients 
are even. In the following result, pi = ini Mr)R+ Pii^) , Po = infMnR+ ^2(2;) and 
nil = mi(M n ]R+), ma = m2(M n IR+). 

Theorem 3 Under the conditions of theorem\^ if moreover a, 6, c, d, e, / are 
even functions and ^ M , then there exists an odd non-trivial solution u e 
i/i(M) X ifi(M) of equation such that r < ||u|| < R. 

Proof: The proof mimics the steps of the proof of Theorem [21 working now 
with the Green's functions for the problem on the half-line M+, with boundary 
conditions Ui{0) = = Ui{+oo). The operator T : BG(R+) x BC{R+) 
BG{R+) X Bc{r+) is defined by 
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Ti{u){x)= Gi{x,s) {b{s)u2{s) + c{s)F{ui,U2)ui{s))ds (19a) 
Jr+ 

T2(u)(x)=/ G2ix,s){e{s)ui{s) + f{s)H{ui,U2)u2is))ds (19b) 

Jr+ 

The adequate cone is 

P = {u G BC{R+) X BC{R+) : Mi(0) = ^2(0) = 0, 

Ui(x),U2(x) > 0,Va; G M+, min Ui(x) > mj»nll'"ill r • (20) 

xeMnm+ ^ un n j 

Then T{P) C P and T is a continuous and completely continuous operator, 
since Proposition [1] can be applied to functions u defined only in M+ and 
such that u{0) = by simply extending as the zero constant function on the 
negative axis. Finally, the sets Qi and ^2 are defined again as the open balls 
of radius r and R, respectively. Everything works in the same way so the 
repetitive details are omitted. In conclusion, we obtain a positive non-trivial 
solution u e H\R+) x H^{R+) such that mi(0) = ^2(0) = and the odd 
extension gives the desired solution.a 

Note that the assumptions ^ M is necessary in order to have rrii 7^ 0, 
1712 7^ 0, which is a key point in the proper definition of the cone. 

Finally, let us note that an analogous of Corollary [T] also holds for this case. 



5 Branches of solutions in situations of physical interest 

As it was said in the introduction, there are two situations arising in physical 
applications described by Eqs. (l5al) and (l5bl) . Our results can also provide 
some information on the localization of the solutions that can be of interest 
in the study of branches of solutions in systems controlled by parameters. As 
a basic example, we consider the system (jll) with the nonlinear contribution 
given by Eq. fl5a|) . i.e. 

— u"{x) + a{x)ui{x) = h{x)u2{x) + \c{x)u\{x) (21a) 
-u'2{x) + d{x)u2{x) = e{x)ui{x) + Xf{x)ul{x) (21b) 

where A > 0. 

Corollary 2 Let us assume the conditions of Theorem\^ Then, for all X > 
there exists a positive solution ux = (mia,W2a) £ H^{M) x H^{R) of l{21a\) . 
Moreover 

lim ||u\|| = +00, lim ||ua|| = 0. (22) 
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// moreover a, b, c, d, e, f are even functions and ^ M, there exists a second 
branch of odd solutions Ua € H^{M.) with the same limiting properties. 

Proof: The application of Theorem [2] requires the existence of r^, R\ such 
that 



1 — max ( / Gi{x, s)b{s)ds, / G2{x, s)e{s)ds 



M 



A max 



Gi{x, s)c{s)ds, G2{x, s)f{s)ds 



M 



1 -1 



AmoPo max / Gi(x, s)c(s)ds 
x&M Jm 



(23) 



Since Jj^j Gi{x, s)b{s)ds, f^j G2{x, s)e{s)ds < 1 for all x G M, such r^, R\ exist 
and can be chosen so that 



lim Tx = +00, lim Rx = 0. 



(24) 



Hence, we obtain a branch of solutions such that rx ^ ||ua|| ^ Rx, and 
now a passing to the hmit finishes the proof. The arguments for the branch of 
odd solutions are analogous. ■ 

Thus, we obtain a bifurcation from infinity when A — 0^. 

We can also include the linear part under the effect of the parameter A, i.e. 



— u'l{x) + a{x)ui{x) = A (b{x)u2{x) + c{x)ul{x) 
—U2{x) + d{x)u2{x) = A (^e{x)ui{x) + f{x)ul{x) 

where, now, A G (0, m) is a positive parameter, with 



m 



max 



Gi{x, s)b{s)ds, / G2{x, s)e{s)ds 
Jm 



1 -1 



(25a) 
(25b) 



(26) 



xeM \Jm 

Then, we obtain the following result. 

Corollary 3 Let us assume the conditions of TheoremlM Then, for all A G 
(0,m) there exists a positive solution ux = {uix,U2x) G H^(M) x H^(R) of 
^25 a\) . Moreover 

lim IIuaII = +00 (27) 

A^0+ 

// moreover a, b, c, d, e, / are even functions and ^ M, there exists a second 
branch of odd solutions Ux G H^{M.) with the same limiting properties. 

By using the same argument as in Corollary [H one get an explicit lower bound 
for m. 
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